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Abstract: We study the time-dependent solitonic gauge fields in scalar QED, in which a 
charged particle has the energy of reflectionless Posch- Teller potential with natural quan- I ^ 

turn numbers. Solving the quantum master equation for quadratic correlation functions, we 
find the exact pair-production rates as polynomials of inverse square of hyperbolic cosine, ^ ^ 

which exhibit solitonic characteristics of a finite total pair production per unit volume and 
a non-oscillatory behavior for the entire period, and an exponentially decaying factor in 
asymptotic regions. It is shown that the solitonic gauge fields are the simplest solutions of 
the quantum master equation and that the back-reaction of the produced pairs does not 
destabilize the solitonic gauge fields. 
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1. Introduction 



Schwinger showed that a constant electric field could produce pairs of charged fermions in 
quantum electrodynamics (QED) and of charged bosons in scalar QED Parker showed 
that an expanding spacetime could also produce particle pairs with opposite quantum 
numbers Under background gauge fields or curved spacetimes the out-vacuum can 
differ from the in- vacuum and the nontrivial Bogoliugov transformation between two vacua 
leads to pairs of the out- vacuum particles in the in- vacuum or vice versa ^ . The recent 
proposals to directly test Schwinger pair production and vacuum polarization effects in 
ultra-strong laser sources js], ^, ^ § and possible astrophysical applications have 
boosted studies on vacuum instability and pair production |11, [l^, |l^ 



14, 15, 16, 17, 18 



19, 20, 21, 22]. Several field theoretical methods have been elaborated to calculate the 

m,m (for 



pair-production rate in various electric field backgrounds [23, 24, 25, |2^, 27, |2^^ 

The author and Schubert have recently introduced 
32| , which relates the quantum master equation for 



review and references, see refs. |31, |T^^ 
a time-dependent solitonic gauge field 
quadratic correlation functions to the KdV equation, whose soliton solution determines 
the pair-production rate and the other correlation functions. One interesting feature of the 
solitonic gauge field is that the pair-production rate is determined by a single hyperbolic 



function for the entire period, which exponentially increases from and decreases back to 
zero in a symmetric way without any oscillatory behavior throughout the evolution. 

In this paper we find a more general class of the so-called "solitonic gauge fields" 
that lead to the pair-production rate and other correlation functions having the "solitonic 
characteristics" not in the sense of traveling solitons for the KdV equation but in sense 
of finiteness of total number of pairs produced for the entire period per unit volume, 
an exponentially decaying behavior in asymptotic regions and a non-oscillatory tailing 
behavior. Further the back-reaction of produced pairs should not destabilize the solitonic 
gauge fields. These solitonic gauge fields have the form of Posch- Teller potential in quantum 
mechanics, which also provides the reflectionless potentials whose bound solutions lead to 
solitons for the KdV equation |^3| . A massive scalar field in the global geometry of de Sitter 
spaces has the Posch- Teller potential for each spherical harmonics, in which particles are 
produced only in even dimensions but not in odd dimensions There the 

pair-production rate, which is obtained by comparing the out-vacuum in the future infinity 
with the in-vacuum in the past infinity, is an asymptotic result and does not give any 
kinematic information on particle production during the evolution. It is argued that pair 
production in de Sitter spaces may cause the decay of the vacuum energy |3^ . 

We show that a charged boson with the Posch- Teller potential for the squared energy, 

with the free energy ujq, has the solitonic pair-production rate for a natural number p, which 
is a consequence of the reflectionless scattering. The gauge field strength is determined by 
the parameter p. The massive scalar field in the global covering of de Sitter spaces has the 
correspondence of the parameters 

I'^i^^Y-ciy, ,^,+r>=^+i(^-l) (1.2) 



where H is the Hubble constant, d is the spatial dimensions, and —A is the eigenvalue 
the Laplace operator on S"^. The main difference of de Sitter space from QED is that 
the mass of the particle is the free parameter while all the other parameters, such as the 
Hubble constant, the spatial dimensions and the eigenvalues of the Laplace operator, are 
predetermined. As the parameter p can be controlled as the interaction strength in QED, 
it is legitimate to call the corresponding fields as the "solitonic gauge fields," though all 
the results of this paper hold directly for de Sitter spaces. 

The method to be employed for finding the pair-production rate is quantum kinemat- 
ics for the correlation functions of the spectrum generating algebra of the annihilation and 
creation operators for the charged bosons in a gauge field background. The Klein-Gordon 



equation in the gauge field has the SU{1, 1) algebra [32|. This means that the correlation 
functions carry the same information as the mode solutions or the propagator of the Klein- 
Gordon equation. Though the pair-production rate may be found from the mode solutions, 
which is the case of solitonic gauge fields, the advantage of quantum kinematic approach 
is that the production-rate and the correlation functions for quadratic annihilation and 



creation operators obey a system of linear differential equations with time-dependent co- 
efficients, which can be integrated from the initial moment in which no particle is present 
and the gauge field vanishes asymptotically. 

Solving the quantum master equation, we find the exact pair-production rate in the 
background ( |f .f| ), which exhibits the solitonic characteristics. The solitonic solutions may 
shed light on two important questions in particle production: the local nature of and the 
oscillatory behavior of pair production. The effective action carries both the global and 
the local information of the background fields either through the eigenvalues of the field 
equation or through the real-time propagator |4^] or through the Bogoliubov coefficients 
from the asymptotic mode solution [41, 42, 43]. In the in-out formalism the vacuum per- 
sistence (twice of the imaginary part of the effective action) relates to the pair-production 
rate at one-loop |^, Q 

2Im(£S(t)) = J]ln(l+AAK(t)), (1.3) 

K 

where K denotes all quantum numbers. Therefore a physically interesting question may be 
raised where and when to observe produced pairs. Another surprising aspect is that though 
the Schwinger formula in a constant electric field is constant [Q], the pair-production rate 
in general time-dependent electric fields shows the oscillatory tailing behavior as explained 
in detailed below. 

We also advance an algebraic method to directly solve the quantum master equation. 
The dependent variables for the quantum master equation are the quadratic correlation 
functions of annihilation and creation operators, which are nonlinear functionals of the 
complex mode solution, the first derivative, and their complex conjugates. Therefore the 
quantum master linear equation, equivalent to the first order master equation for three 
independent correlation functions, is a highly nonlinear functional of the mode solutions 
||3^ . However, in terms of the new variable the master linear equation takes a relatively 
simple form with the coefficients of the gauge field and its derivative, which allows an 
algebraic approach to the quantum kinematics. Under the condition that the square of 
the first derivative of the profile function for the gauge field is a polynomial of the profile 
function itself, we are able to find the recurrence relation for the power series solution 
for the quantum master equation. This algebraic method does not assume the profile of 
the gauge field a priori but determines the profile function from the requirement that the 
recurrence relation should terminate for soliton solutions. Remarkably we find that the 
Posch- Teller potential (1.1) with natural numbers p are the simplest solitonic gauge fields 
from the algebraic method. However, the algebraic method leaves open for other solitonic 
gauge fields. 

The solitonic gauge fields (natural number p = n) strongly contrast with other non- 
solitonic gauge fields (non-natural number p) of the same potential shape which predict the 
oscillatory tailing behavior. The oscillatory behavior may be understood from the quantum 
master equation in the non-adiabatic approach [^] 

|(l + 2AAk(t)) = "^^'^"^.^^'"^ fdi\{u:l{i)-u:l{t,m + 2M^{i)) 



X COS 



Wk(io) 



(1.4) 



or the quantum Vlasov equation from the adiabatic approach |44, 45, 47, 48, 12] 



k(0)cos(2 / d^'a;^(C' 

to 



(1.5) 



where ujk{t) is the energy in the time-dependent electric field and uJk{to) denotes the 
initial value. In the asymptotic regions {t = ±oo), eq. (|l.4|) or (^) weights the en- 
velop of exponentially decreasing function with the sinusoidal function dominated either 
by cos{2uji(_{to){t — ^)) in the non-adiabatic approach or by cos{2uji(_{to)t) in the adiabatic 
approach. Therefore the non-oscillatory tailing behavior of solitonic gauge fields is an 
astonishing result since the oscillations should exactly cancel. 

We investigate the back-reaction of the produced pairs by computing the induced 
current and the induced electric field in a one-soliton gauge field. The renormalized current 
and thereby the induced electric field has the leading term proportional to the solitonic 
gauge field background by a factor of the fine structure constant. This means that the 
back-reaction through the induced current and the electric field may be neglected for the 
solitonic gauge field. The full scenario for the back-reaction problem is the quantum master 
equation in the gauge field background plus the induced gauge field and the equation for 
the current, which makes the system coupled and highly nonlinear. We leave it as an 
open question whether the solitonic gauge fields could keep the same shape or whether any 
solitonic gauge field could exist or not. 

The organization of this paper is as follows. In sec. ^ we exploit the properties of 
solitonic and non-solitonic gauge fields. In sec. ^ we find the pair-production rate from 
the mode solutions. In sec. Q we find the power series solutions of the quantum master 
equation for the solitonic and non-solitonic gauge fields. In sec. ^ we solve the quantum 
master equation by an algebraic method and show that the solitonic gauge fields are the 
simplest one allowing polynomial solutions. In sec. ^ we investigate the pair-production 
rate and correlation functions in solitonic gauge fields. In sec. |^we compare pair production 
in non-solitonic gauge fields with that in solitonic gauge fields. In sec. ^ we study the back- 
reaction of produced pairs by computing the induced current and electric field and show 
that the back-reaction is always order of the fine structure constant and thus negligible. 



2. Solitonic and Non-Solitonic Gauge Fields 

The recently introduced time-dependent solitonic gauge field |32] relates the quantum mas- 
ter equation to the KdV equation, whose solution determines the pair-production rate and 
the correlation functions of pair annihilation and creation. We first generalize it to (i) the 
solitonic gauge fields 

qA^^^^^it) = ^^-(kl + ^^^^±^f\ (n = l,2,3,---), (2.1) 



where q and m are the charge and mass and wq = Vk^ + m? is the free particle energy 
in the asymptotic hmits t = ±00. The case of n = 1 is related with the KdV equation 
|32]. For the zero-longitudinal momentum the gauge fields simply become = 
— Y^n(n + l)a;o/('? cosh(a;ot)) and the potential difference is zero between t = —00 and 
t = 00, and the electric field 

Vn(n + lVg sinh(a;ot) 

^{n)\\[t) = TIT — 7: y^-^i 

^ q cosh^(wo*) 

changes the polarity and exponentially decreases in the asymptotic limits t = ±00. The 
Fourier-component of a charged scalar in the solitonic gauge field obeys the mode equation 
[in units of h = c = 1] 

As the solution depends on the free particle energy, we rescale the time as t := oj^t and 
measure the time in the Compton time with the equivalent mass of relativistic energy or 
equivalently set ojq = 1 for simplicity, delete the momentum k unless necessary, and denote 



'^(n)W:=l + ^^^, (n = l,2,3,---). (2.4) 
The next generalization is (ii) the two parameter-dependent gauge fields 

^L)(0:=^' + ^^^, (p/1,2,3,-..) (2.5) 



where / is positive but not necessarily natural numbers. In quantum mechanics the poten- 
tials of the form (2.4) and (2.5) are known as Posch- Teller potential. These potentials have 



a rich structure and a variety of interpretations for the Klein-Gordon equation, as will be 
shown below. 

The energy ( |2.4D or ( |2.5D with a natural number p = n is known as the reflectionless 
potential for the mode equation ( p.3[ ), whose bound solutions for another natural number 
I generate the solitons for the KdV equation [^3|. The reflectionless scattering may be 
understood semiclassically by computing the complex Hamilton-Jacobi action in the phase- 
integral method 0, 21|. Taking into account the four complex turning points, the 



leading term of the pair-production rate is 



M « 4cos2(7rv^n(n + l))e-2'^, (2.6) 



which destructively interferes approximately for large n such that n{n + 1) ~ (2n + l)/2 
or completely interferes away by including the Maslov index contribution 7r/4 to the action. 



The exact solution for the gauge fields (2^) is given by the hypergeometric function |51| 

-ilt 

(^(p;)(t) = ^(1 + e^'Y^'F{p +1 -il^p +1-1 -il- -e2*), (2.7) 
V 2t 



which has the proper asymptotic behavior at t = — oo for the Gaussian vacuum with the 
energy l{ujo): 



-lit 



In the other asymptotic Hmit t = oo, the mode solution takes another form 



(2. 



-ilt 



T{l-il)T{-il) 



21 T{p+l-il)ri-p-il) 



[1 + e-^y+^Fip +l,p+l + il;l+il; -e 



2t 



Jit 



+ 



T{l-il)T{il) 
2ir{p+l)T(-p) 



(1 + e-2*)P+iF(p + l-il,p+l;l-il; -e' 



2t\ 



) 

2.9) 



The first case (i) is when p is a natural number n, for which the negative frequency vanishes 
due to r(— n) = (—1)" x oo. The solution describes a reflectionless scattering over the 
potential, so pairs are not produced in the asymptotic region t = oo. In the second case 



(ii) of a non-natural number p, the mode solution (2^) has both the positive and the 
negative frequencies in eq. (p.9|), implying pair production in the limit t = oo. The mode 



solution may be further written in terms of the associated Legendre function [51, 52] 



¥'(p/)(0 



T{l-il) 
~72~ 



P;'(-tanht) 



(2.10) 



The energy (2.4) or ( p.5| ) has a finite number of bound solutions for positive imaginary 
I = in. The mode solution (|2.7D exponentially decreases at t = — oo while its analytical 
continuation ( |2.9D to f = oo has the exponentially increasing branch, the first term, and 
the exponentially decreasing branch, the second term. The exponentially increasing branch 
vanishes when p — k is a non-negative integer: the number of bound solutions is [p], the 
largest integer no greater than p. For the natural number p = n, the bound solutions are 
the associated Legendre polynomials P^(— tanht) for / = 1,2,--- ,n, which generate the 



n-solitons for the KdV equation [33| 



3. Pair Production in Terms of Mode Solution 



In this section we use the mode solution (|2.7| ) and the other asymptotic form ( |2.9D to 
calculate the pair-production rate in the limit t = oo. However, one may still ask whether 
the pair production is absolutely zero or decreases to zero in that asymptotic limit and ask 
the kinematical information of pair production during interaction, which will be pursued 
in this paper. 

To find the exact quantum states, we may employ the Lewis-Riesenfeld invariant op- 



erators [53] in the Schrodinger picture, which provide the time-dependent creation and 



annihilation operators for particles and antiparticles p4|, B5L 56 



«{pOkW = -^[v'(p0k(0^(pZ)k-V5(p/)k(i)'^(p/)k]> 

a(pi)Ut) = ^[¥'(p/)k(*)^(pOk - '^{pOk(*)'^(pOk]' (3-1) 



and for antiparticles 

^(pO-k(*) = -^['/'(p/)k(i)vrfp;)k-9i(p/)k(i)<^(pOk]> 

ft{pO-k(0 = ^KpOk(*)^(p/)k - 'P*pi)kit)(f>lpi)],]- (3-2) 

Here 4'{pi)\i s-^icl vr^p;)^ are the Fourier component of the field and the momentum, all being 
Schrodinger operators, and (/?(p;)k is the complex solution ( |2.7] ). The mode solution ( p^ ) 
satisfies the Wronskian condition from the canonical quantization 

Wr[(/9(pi)k, (^(p;)k] = V'(pOk(*)'/5(p/)k(*) - V(p/)k(*)'/^(pOk(*) = h (3.3) 
which is equivalent to the equal-time commutation relations 

[«(pz)k(t),a[pi)k'(i)] = [^(pO-k(*).^(pO-k'(*)] = '^(k - k'). (3.4) 

Then the Bogoliubov transformation between the past infinity to = ~oo and the present 
time t is given by 

a(p/)k(-oo) = A^(p/)k(-oo,*)a(p/)k(0 + i^(p/)k(-oo,*)^(p;)_k(^)' 
^(pO-k(-°o) = KpOk(-°°' *)^(pz)-k(*) + '^(pOk(-oo> *)a(pOk(i)> (3-5) 

where 

^(pOk(-oo,t) = iWr[99^p^),^(-oo),v?(p;)k(t)], 

'^(p/)k(-oo,t) = iWr[(^^p;)k(-oo)'V2(p0k(*)]- (3-6) 

The Wronskian (|3.3| ) makes the Bogoliubov coefficients satisfy the bosonic relation 

lAt(pOk(-oo,t)|^ - |z^(pi)k(-oo,t)p = 1. (3.7) 

The time-dependent ground state |0(p;)k)i)) a Gaussian state, is nullified by the anni- 
hilation operators for a given k 

a(pOk(^)|0(pi)k,i) = ^(pO-k(*)|0(p/)k)*) = 0) (3.8) 
and the time-dependent vacuum is the tensor product 

\%i).t) =\{\%i)^,t). (3.9) 

k 

The mean number of pairs carried by the time-dependent vacuum as measured by the 
in-vacuum number operator is given by 

(0(p/),i|a|p;)k(-oo)«(pOk(-oo)|0(p,),t) = \y{pi)i,{-oo,t)\^, 
(0(p;),*l&(pi)_k(-oo)V)-k(-oo)|0(pi),t) = |j^(pi)k(-oo,t)|^. (3.10) 

Similarly, the mean number of pairs contained in the in-vacuum now as measured by the 
number operator at t is the same as eq. ( ^.lOl) since V(pi)\i{t.,—co) = '^(pi)k(~°*^' ^^^^ 



equal number of particles and antiparticles are produced due to charge neutrality of the 
in-vacuum and the electric field. Interpreting as the number of pairs the expectation value 
of the time-dependent number operator in the in-vacuum or the number of in-vacuum 
particles in the time-dependent state is a subtle issue (see refs. [47, [l^). In the adiabatic 
case in which |ti;^/a;^|, |ci)/a;^| <^ 1, we may use the adiabatic solution to count the adiabatic 
particle number and the number of pairs ( 3.1C| ) approaches the adiabatic one. The 
gauge potentials ( p.4| ) and ( |2.5[ ) satisfy the adiabaticity, which will be discussed in detail 
in sec. The pair-production rate ( 3.10| ) from the solution ( ^.7| ) is 

1 



= -^lApiM , (3.11) 

where ^(p/) is the time-dependent amplitude 

^w(i) = (1 + e^y+'Fip + l-il,p + l-l-il; -e2*). (3.12) 
iboo, the pair-production rate for the solitonic gauge field is 



In the asymptotic limits t 
approximately given by 



n(n-l-l)\2 1 



21 



1 



4(l + /2)cosh^t' 



(3.13) 



The pair-production rate in the non-solitonic gauge field still has the same form ( |3.13 ) in 
the limit t = — oo, but the other limit t = oo will be separately treated in sec. ^. 

The physical meaning of I in eq. ( |2.5D is the squeezed vacuum of the Minkowski one 
for 1 = 1. The squeezed vacuum is the Gaussian vacuum with the energy /(wq) while the 
Minkowski vacuum has the minimum energy 1{ujo) and the uncertainty 1/2 p9| , |60t| . Com- 
paring the asymptotic solution (|2.8| ) for / 1 with that for / = 1, we find the Bogoliubov 
coefficients 

M(nZ)(n) = i^^[^*n)(i)^V>{nl){t)] ^ ^ ^ 



-i(l-l)t 



l-l 



-i(/+l)t 



(3.14) 



The Bogoliubov transformation ( [3.14] ) implies that the vacuum state for I is the squeezed 
vacuum of the Minkowski one in the asymptotic limit t = -co Ell] 



|0(n/),i) 



exp 



n(a(n)(t)fe(n)(i)e'' -a;„)(i)^(„)(i)e 



~2it 



X exp 



^(^ - 1)^( «(„) Wa(n)(i)) + hn){t)b\^){t) 



|0(n))*)) 



where 



cosh ri 



l + l 



sinh ri 



I - 1 



(3.15) 



(3.16) 



2^^l' 

The probability for the squeezed vacuum to be in the number state | , t) of the Minkowski 



vacuum is 



Pk = \ {k(n),i\0{ni),i)? 



tanh^''' ri 



41 



coshV, (/ + 1)2V/ + 1 



/ - l\2fc 



(3.17) 



Thus the squeezed vacuum has the mean number of pairs with respect to the Minkowski 
vacuum and its Fock space 



J^(nl)(n) = ^ = — = \v(^nl){n)? ■ (3.18) 



Al 

k=0 



4. Series Solution for Master Equation 

The hnear combinations of quadratic annihilation and creation operators for each compo- 
nent of the Klein-Gordon equation, not including interaction terms, in an electromagnetic 



field background or in a curved spacetime, constitute the SU (1,1) algebra ||3^. The evo- 
lution of these generators carry the same information as the mode solutions in sec. |3|. 
The dynamical equation that determines the pair-production rate from an initial parti- 
cle distribution is of the primary concern. For that purpose, we note that the first order 
time-derivative of the pair-production rate 

X{pi)it) ■ = (0{pO'*I^(p/)(~°°)"{pO(~°°) + V)(~°°)^(p/)(~°°)|0{pO'*) 

= l + 2%;)(t) (4.1) 

is related with the other correlation functions 

Z{pi)it) ■■ = (0(p,),t|a(p;)(-oo)6(p/)(-oo) -Fa|p;)(-oo)6{p,)(-oo)|0(p/),t). (4.2) 
The unity in eq. (O) is the zero-point energy of quantized field, which relates to renor- 



malization of the vacuum energy. In fact, the quantum master equation consists of coupled 
linear differential equations 

. _ p{p + 1) 

• _ P{P +1) , f^r . P(P+1) Y^ 

The initial conditions are such that no pair is present and the correlation functions vanish 
in the asymptotic region t = — oo 

X{-oo) = l, y(-oo) = Z(-oo) = 0, 

X{-oo) = y(-oo) = Z(-oo) = 0. (4.4) 

The quantum master equation ( [4. 3]) is equivalent to the linear equation [32] 

flM + 4(/2 + ^^)^ + 2(-^^^)F(,, - = 0, (4.5) 

dt^ V cosh^t/ dt Vdtcosh^t/ ^^'> dt P coshH ^ ^ 



where 

F{pi)it) = \ f diY^^DiO- (4.6) 
The pair-production rate is then given by 

X^it) = 1 + mvDiO^-^^, (4.7) 
and the other correlation functions are given by 

z^it) = - |*^de%,)(e)(2/2 + ^±i)). (4.8) 

Closely inspecting eqs. ( |4.7D and ( |4.8| ) and using another relation |32| 

-^ii) = ^2 Ho(*)^(*) - 2 / ^ dCmHo^^) (4-9) 

leads to the conserved quantity for the solution F to eq. (|4.5D 

X^{t) -Y'^it) - Z^{t) = 1. (4.10) 



The particular form of eqs. ( [4. 3D or ( [4.5[ ) or ( [4.7] ) and ( |4.8| ) suggests the power series 
solution 

FiPi)it)=T.C(pi)k{^^y, (4.11) 

k=l 

and 

-^(ni) jt) = -P(n;)fc ( ^^^-^^2 ^ ) ' -^(n/)0 = 1. (4.12) 
k=0 

The recurrence relation for C(^pi^k can be found, for instance, from eq. ( [4.5| ) 

_ (2fc + l)(p-fc)(p + fc + l) 
^(pOfc+i 2(A; + l)((A: + l)2 + ;2) ^W^' 

_ P(P + 1) , . , ON 



and the coefficients are found from eq. ( [4. 7] ) 



_ p(p+l)fc 

^(P')fc - ^ ^ ^ '^(pOfc- (4.14) 
The coefficient C(p;)^ may be written in a compact form 



-10- 



5. Soliton Solutions by Algebraic Method 



In this section we directly find the soUtion solutions by an algebraic method. We shall 
confine our study to the gauge field leading to the squared energy 

Lo\t)=ioi{l^ + Uou{t)), (5.1) 

where the profile function u{t) has a finite integral over (—00, 00), and u{t) and its derivative 
u{t) approach zero sufficiently rapidly in the asymptotic limits t = ±00. And we further 
assume that u{t) has the following property 

t^2(t) = ^Pfcn^(t). (5.2) 

k=l 

There is no constant term since u{±oo) = n(±oo) = 0. The solitonic gauge field u = 
1/ cosh^(wot) has the property (|5.2|). In the rescaled time t := uj^t, the quantum master 
linear equation ( |4.5| ) reads 

__ + 4(/^ + Uou)- + 2U,-T - = 0. (5.3) 

We search for the polynomial solution of the form 

N 

F(t) = Y,Cku\t), (5.4) 

fc=i 

which shares the same solitonic characteristics as u. In order for the series solution to be a 
polynomial for some finite A^, the negative contribution to from d?F/dt^ should cancel 
the positive contributions from udJ- /dt and uJ-", otherwise the recurrence relation will have 
all positive coefficients and will become an infinite series. For this reason, we assume an 
ansatz for the gauge field 

i? = av? + I3u^ , (a > 0), (5.5) 

which has the soliton solution 

u{t) = " ^ . (5.6) 

;3cosh2(^) 



Furthermore, the master linear equation ( |5.3D together with the ansatz ( |5.5| ) leads to the 
simplest linear recurrence relation 

^ (2fc + l)[;gfc(fc + l) + 4[7o] 
2(fc + l)[a(A: + l)2 + 4/2] 

= /2(a + 4/2)- (5-7) 

The series terminates when Uq/ (3 = —N{N + l)/4 for some natural number N . We recover 
the solitonic gauge fields in sec. |2| by choosing a = 4 and (3 = —4, and the recurrence 



relation ( |5.7D reduces to ( f4.13| ) for the natural number p = N for soliton solutions in sec. 
i 

On the other hand, another simple ansatz, for instance, -u^ = au + fSu^, leads to a 
sinusoidal gauge potential, not a solitonic gauge field. Still another ansatz 



leads to a rational soliton solution |33] 



(/3<0), (5.9) 



(at)2 - 4/3 

but the recurrence relation involves three terms 

(2fc + l)[afc(fc + l) + 4[/o] ^ 13 
= 8P{k + l) 

Ci = % (5.10) 

The finite series requires /3 = and C/q/q = —N(N + l)/4 for some natural number A^, 
which makes the solution (^^) ill behaved near t = 0. Solving (|5l0| ) and summing the 



series ( ^.41 ) to get a closed form is beyond the scope of this paper, which will be left as an 
open problem. 

6. Pair Production in Solitonic Gauge Fields 



In this section we consider (i) the solitonic gauge fields (|2.4| ) , where p = n is a natural num- 
ber. The solitonic gauge field corresponds to the refiectionless scattering of the incoming 
positive frequency solution and predicts no pair production in the asymptotic limit t = oo. 



The series ( 4.11 ) becomes a polynomial solution 

k=l 

since C(^ni)k — foi' k > n + 1. Note that the time is rescaled as t := loqI. The coefficients 
change signs alternatively and the solution exponentially decreases in the asymptotic re- 
gions (t = ±00), exhibiting solitonic characteristics. In fact, the solitonic solution is a finite 
expression for the entire period for another exact formula ( ^.11[) from the mode solution, 
which is an infinite series. Though eq. ( 3.11| ) is equal to eq. ( |6.lD , any transformation 



formula of the hypergeometric function is not available to our knowledge 52 1. The 
solitonic solutions provide the global kinematics for the pair-production rate and show the 
temporal nature of pair production. The first four solutions of -^"(71/) for / = 1 are given by 

^"■'<" = ^i?7- <«-^* 

3 9 9 

^ 5 27 63 441 

^^^^ 2 cosh^ t 4cosh^t ^ 8 cosh*' t 136 cosh^t' 




Figure 1: The functions -^(711) and the pair-production rates Af(ni) are plotted against t{:— ujQt): 
the solid line for the dotted line for i^(2i), the dashed line for i^(3i), and the dash-dotted line 

for i^(4i) [left panel]; the solid line for A/'(ii), the dotted line for A/(2i), the dashed line for A/(3i), 
and the dash-dotted line for A/'(4i) [right panel]. 



The one-soliton solution ( |6.2D is the same as eq. (25) of ref. [32[. The profiles are relatively 
simple and their variations depend on quantum number n, as shown in fig. |l|. 
From eq. ( [4.7D the pair-production rate is given by 



n(n + l)^ k ^ / 1 



fc=i 



The first term is the leading term ( 3.13 ) from the mode solution in the asymptotic regions 



{t = ±00). The pair- production rate in the first four solitonic gauge fields for / = 1 are 

AA(3l^(^) = V + (6-9) 

2cosh^t cosh^t lecosh^t ^ ^ 

. , , , 25 45 945 441 

J\f(4i)(t) = J ^ + ^ T^-. (6.10) 

^ ' 2cosh^t cosh^t lecosh^t ITcosh^^t ^ ^ 

The pair-production rates in fig. || exhibit the solitonic characteristics. They increase as 
quantum number n, the strength of interaction. The pair-production rate from the master 
equation exponentially decreases to zero, does not have any oscillatory tailing behavior, and 
exhibits the solitonic characteristics. The momentum distribution of the pair-production 
rate has a substructure for time less than one Compton time t{:= mt) = 1, as shown in 
fig. |2|. And the substructure appears within the domain of a few Compton times and 
wavelengths larger than a few Compton wavelength, as shown in fig. ^. In fact, the 
substructure depends on the combination of ujQt = ^Jl + {k/mY{mt). 



We now turn to the second case of (2^) with natural numbers n and / > 2. The mode 



solution ( |2.7D does not have the negative frequency branch even for any / at i = 00 as 
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Figure 2: The momentum distributions A/'(„i)(fc) are plotted against k{:~ k/m): at t{:— nit) — 0.1, 
the solid line for A/'(ii)(fc), the dotted line for A/'(2i)(fc), the dashed line for A/'(3i)(fc), and the dash- 
dotted line for A/'(4i)(fc) [left panel]; at t{:= mt) = 1, the solid line for the dotted line for 
A/'(2i), the dashed line for A/'(3i), and the dash-dotted line for A/(4i) [right panel]. 





Figure 3: The time-momentum distributions J\f(2i){t,k) are plotted against t{:= mt) and fc(: 
k/m) [left panel] and A/'(4i) (t, fc) [right panel]. 



shown in sec. |2|. The master linear equation ( [4 .51 ) still has the polynomial solutions 
with (4.15). In fig. ^ we compare -^"(44) with -^"(42) and also compare -^"(32) with -^(33^. The 



non-natural number / drastically changes the profiles of -F(4) and -F(3) with / = 1 and other 
natural number I. 



7. Comparison with Non-Solitonic Gauge Fields 

Finally we consider the mode equation in the non-solitonic gauge fields with positive non- 
natural number p 



(7.1) 
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Figure 4: The functions Ft^nX) are plotted against t{:— uot): the sohd hne for i^(44) and the dotted 
hne for ^"(4|) [left panel]; and the solid line for i^(32) and the dotted line for Ff^^^-^ [right panel]. 



One is tempted to use the series solution (4.11), which becomes an infinite power series 
since the coefficients do not terminate for non-natural number p. The series solution ( 4.11| ) 
can be written as 



[p] 

k=l 



1 



cosh^ t 



+ 



E 

k=[p]+l 



{pl)k 



cosh^ t 



(7.2) 



where the first sum is finite and alternating while the second sum is non-alternating and 
infinite. The series (^^) converges strongly near the limit t = — oo and gives an accurate 
result. The first term is the leading term for the pair-production rate (|3.13| ) from the mode 
solution. As the coefficients ( f4.15 ) do not change the sign and approximately approach a 
constant value Cn times 1/k for k > N for some large N, the series solution for large k 
converges asymptotically to a logarithm 



F(p,)(t)«-C;vln 1 



1 



cosh Vcosh t 



1 



N 



(7.3) 



Therefore the series solution cannot be relied on when coshf ~ 1 and near t = 0. This 
implies that the formula ( [4.7D for pair-production rate may not be used beyond the limit in 
which the series ( [7.2[ ) converges since it is not justified to evolve ( [7. 2D from one asymptotic 
limit t = — oo to the other limit t = oo through t = 0. Thus the series solution to the master 
equation does not provide useful information on the global kinematics for pair production 
beyond t = more than the mode solution, as shown below. 

The exact solution (p.7|), which is valid for the entire period, can be transformed into 
another form useful in the other asymptotic limit t = oo 



-ilt 



V2 

Alt 
. ^ i- 

H — ^e*2 
\/2 



sin^(7rp) cosh^(7r/) -|- cos^(7rp) sinh^(7r/) \ 1/2 



sinh (vr/) 



B{pi){t) 



sin(7rp) 
sinh(7r/) 



(7.4) 



where 'dr is the phase of the gamma functions of (|2.9| ) and B(,pi-^[t) is another time- 
dependent amphtude 



B(pi){t) = (1 + e-2*)P+iF(p +l,p + l + il;l + il; -e'^'). (7.5) 
A straightforward calculation gives the first two leading terms of the pair-production rate 

_ ( sin(7rp) \ ^ P(P+ 1) -2A n P(P+ 1) / sin(7rp) 

■^^pOW - ^sin^vrO/r 1 + /2 ^ J ^(l + /2)i/2Uinh(vr/)y' 



sin2(7rp) cosh^(7r/) + cos2(7rp) sinh^(7r/) \ 1/2 



sinh2(7r/) 



j sin(/t)sin(/t+ ^ -i?r + i?/)e~^*, 

(7.6) 



where 



The frequency of 21 in the oscillatory factor of the second term is predicted by the non- 
adiabatic and the adiabatic quantum master equations ( |l.4D and (^]^). For a natural 
number p = n, the rate i \l.6\) vanishes as expected. 
For the squared energy of the form 



u^\t) = Vo'+ " (7.8) 
cosh [oJot) 

where Uq is an interaction strength, we identify the parameters 



Noting again t := oj^t and ojq = ^/m? + k2, the constant term in eq. ( |7.6| ) depends on the 
momentum through the sinusoidal function in a nontrivial way. Fig. |5| shows the momen- 
tum distribution of the pair-production rate, which has a substructure for small momentum 
but exponentially suppressed for high momentum. Furthermore, the decaying terms also 
depend on the momentum both through the amplitudes and through the time-dependence. 
The rich substructure of momentum distribution distinguishes the non-solitonic gauge fields 
from the solitonic gauge fields. 

8. Back-Reaction Problem 

In this section we investigate the back-reaction of produced pairs on the solitonic fields. 
The current of a charged scalar particle in an electromagnetic field is the expectation value 
with respect to the vacuum with quantum numbers n and / 
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Figure 5: The momentum distributions oi Af{k) = (sin ttjj/ sinliTrZ)^ for Vo/m = 1 against k{:= 
k/m): the solid line for Uo/rri^ = 1, the dotted line for Uo/m? ~ 4, the dashed line for Uo/m? — 9, 
and the dash-dotted line for XJ^jm} = 16 in the range of — 3 < fc < 3 [left panel] and in the range 
of 3 < A; < 30 [right panel] . 



The first term is the vacuum energy of the quantum field, which has to be regulated away 
through renormalization, and the pair-production rate M{t) and the correlation function 
Z{t) of eq. ( [4. 81) are finite for solitonic gauge fields as shown in sec. To understand 
the time-dependent of the induced current, we consider the simplest case k = for which 
ujQ = m and t := mt. The induced current (^) leads to a linear equation for the induced 
electric field 



EUt) = -{Ji\{t)) 



2q^n{n -h 1) 1 /I 



I cosh t\2 

The induced current, for instance, for one-soliton (n = 1) gauge field for arbitrary I is 



•2) 



£^in(t) 



1 



1 



2(1 + /2) 

1 



1 + 



sinh t 



+ 



1 



sinht 



1 + /2 16/2(1 + /2; 



16/2 7cOsh2t 16/2(1 + /2)^.oghn 

arctan(sinht) 



(8.3) 



As arctan(sinht) = ib7r/2 at t = ±00, we regulate away the last term through renormal- 
ization of the vacuum energy and the charge. Then the ratio of the renormalized induced 
current to the solitonic field background is 



TTiren 



1 + 



+ 



1 



1 



4) 



Sex(t) /(1 + /2)V^ ' 16/2 ' 8/2cosh2t 

and for the electron-positron pairs is order of the fine structure constant (e^) in the units 
oi c = h= 1. Therefore the back-reaction is small compared with the background solitonic 
gauge field. Including the momentum integration, the rate of the renormalized induced 
electric field for one-soliton gauge field is 

d^k 



2g 

/(l + /2)y (27r) 




n(n -|- 1) / 



1 



+ 



cosh (cjot) ^cosh (woi) 4/2 cosh (oj^t) 



5) 



where the dots denote a subtraction scheme for renormahzation. We argue that the back- 
reaction of produced pairs is smaller by the fine structure constant than the solitonic gauge 
field background and can be safely neglected. It remains open to show the existence of 
solitonic gauge fields including the back-reaction, which is the beyond the scope of this 
paper. 



9. Discussion 



One may ask the physical meaning of the pair-production rate during the interaction with 
the electric field. It is defined either by the number of pairs carried by the time-dependent 
vacuum as measured by the in-vacuum number operator or by the number of pairs con- 
tained in the in-vacuum as measured by the time-dependent number operator, both of 
which give the identical result because of the Bogoliubov transformation and the inverse 
transformation ( |3.5| ). In the in-out formalism, the pair-production rate is determined by 
the Wronskian ( |3.6D between the mode solution for the in-vacuum and the mode solution 
at arbitrary time. Does the pair-production rate refer to real pairs to be measured by 
a detector? In the solitonic gauge fields in which the pair-production rate exponentially 
increases from and then decreases to zero in a symmetric way, do the pairs produced during 
the evolution really annihilate themselves in the end? If not, where do the pairs disap- 
pear? This question is rooted on the fundamental question for the particle concept when 
its quantum field depends on time in a nontrivial way, which has been frequently raised 
for expanding spacetimes. 

It is not the purpose of this paper to answer the fundamental question of what is the 
particle concept when its quantum field varies in time. Instead, we pursue a pragmatic 
point of view by searching for regions in which the particle concept and particle number 
makes a sense, at least in the adiabatic sense. In the regions in which the adiabaticity 
holds, we may use the adiabatic mode functions and the adiabatic particle number ||5^. In 
the case of solitonic gauge fields the adiabaticity requires that 

^{nl) _ {n{n + l)sinht)^ 



wf^,) (/2cosh2t + n(n + l))S 
^{nl) _ n{n + 1) cosh^ t 



iof^i^ (/2 cosh^ i + n(n + 1))^ 



n(n + 1)1 ^ +Z^(2cosh^t-3) (9.1) 

cosh t ) 



should be smaller than unity. This condition is always satisfied except for the region 
\t\ < 4 for small natural number n. Even in this limited region the ratios are still less 
than 1/2 and the adiabaticity holds in a broad sense, as shown in fig. ^ In particular, 
the ratios are exponentially suppressed in the asymptotical regions |t| ^ 1 and particles 
are well-defined there. The larger quantum number n corresponds to higher energy than 
the smaller quantum number and thus suppresses further the ratios. The pair-production 
rate begins from the region t = — oo in which the particle number is well-defined, passes a 
region around t = in which the adiabaticity roughly holds, and finally settles down in the 
region t = oo. Therefore we may interpret the pair-production rate as the adiabatic particle 
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Figure 6: The adiabaticity parameter R{t) :— (^)(„;)/i^('„i) is plotted against t{:= cuot): the sohd 
hne for n ~ 1, the dotted hne for n — 2, the dashed hne for n = 3, and the dash-dotted hne for 
n 4 [left panel]. And another adiabaticity parameter S{t) :— ujni / to^j^i-^ is plotted: the solid line 
for n ~ 1, the dotted line for n = 2, the dashed line for n = 3, and the dash-dotted line for n — 4 
[right panel]. 

pairs or in the literal sense as the exact correlation function for the number operator that 
extrapolates the two asymptotic regions. 

It now remains the question of how the pairs produced during the interaction behave 
and disappear such that the pair-production rate vanishes in the other asymptotic limit. In 
the solitonic gauge fields pairs are effectively produced roughly within the interval |t| < 4, 
i.e., eight Compton times (see fig. |l]). During the first four Compton times particles 
accelerate to the positive direction and antiparticles to the negative direction whereas 
during the next four Compton times the directions of particles and antiparticles reverse 
because of the polarization change of the electric field. The pairs that are produced during 
1 < \t\ < 4 may be regarded as real entities since they have enough time to separate over 
one Compton wavelength. However, the pairs produced during jtj < 1 are likely to cancel 
each other due to the polarity change, in a similar way as the virtual pair production and 
subsequent annihilation. Without any detection, the particles keep quantum coherence 
with antiparticles and annihilate during the opposite polarity of electric field. However, 
any detector to measure the pairs may break the quantum coherence and the pairs may 
not completely cancel each other. To install a detector, one has to include an interaction 
between the detector and the charged field, which makes an interaction theory and goes 
beyond the scope of this paper. 

Finally, we compare the solitonic gauge fields with other temporally localized electric 
fields. The frequently studied model is the Sauter electric field 




Eo 





with the gauge field 




(9.3) 



The electric field acts effectively during the period \ t\ < 3r and does not change the polarity. 
Hence a charged particle gains the potential difference of 2qEQT between two asymptotic 



regions and the out-vacuum differs from the in- vacuum, implying pair production |2^, 40 1. 
In fact, the Sauter gauge field ( p^ ) does not satisfy the condition A||(oo) = for solitonic 
gauge fields in sec. ^ which is needed for finiteness of pairs produced for the entire period 
and other characteristics. Another localized electric field 



and its gauge field 



shares some property with the solitonic gauge field. The squared energy 

u,\t)=(ku+^Y + kl+m' (9.6) 



cosh t / 

in terms of the profile function u = 1/ cosht can be transformed into the relation (5.2), but 
the recurrence relation involve more than three terms, which does not necessarily guarantee 
solitonic solutions, except for the zero longitudinal momentum. 

10. Conclusion 

We have studied Schwinger pair production in gauge fields that lead to the energy of Posch- 



Teller potential (1.1) for a charged boson. The gauge fields are classified into two categories: 
the solitonic gauge fields with natural numbers and the non-solitonic gauge fields with non- 
natural numbers. The pair-production rate from the mode solution in a solitonic gauge 
field exponentially decays in the asymptotic regions while a non-solitonic gauge field leads 
to a constant pair-production rate in addition to exponentially decaying terms in the future 
infinity. Solving the quantum master equation, we find the exact pair-production rate in 
the solitonic gauge field which is a polynomial of inverse square of hyperbolic cosine and 
exponentially decays in a symmetric way in the asymptotic regions. Further, the pair- 
production rate in the solitonic gauge field exhibits a non-oscillatory tailing behavior and 
a simple momentum distribution through the entire period. None of these characteristics 
is shared by the non-solitonic gauge fields. The exact information of pair production in 
solitonic gauge fields may shed light on understanding how the vacuum structure changes 
and the produced pairs behave during the interaction with electromagnetic fields. 

Schwinger pair production in time-dependent electric fields can be obtained either by 
finding the Bogoliubov coefficients from the mode equation or by using the real-time prop- 
agator. In this paper we have used the quantum master equation for the three correlation 
functions of quadratic annihilation and creations operators: the number of pairs and the 
correlations for pair production and pair annihilation. The quantum master equation al- 
lows us to directly integrate the pair-production rate from gauge fields only and is easy 
to implement the condition for solitonic characteristics. For that purpose, first we have 



extended the one-solitonic gauge field in ref. [32| to multi-solitonic gauge fields of the same 
potential shape and then solved the quantum master equation to find the pair-production 
rate with solitonic characteristics. Second, we have solved the quantum master equation by 
an algebraic method and showed that the solitonic gauge fields (1.1) with natural number 
p are the simplest ones that satisfy the property imposed on the gauge fields. The algebraic 
method may open a possibility for other solitonic gauge fields. 

In the gauge fields (1.1) the pair-production rate has the meaning of adiabatic particle- 
antiparticle pairs produced during the interaction of the electric field since the mode equa- 
tion satisfies the adiabaticity for the entire period as shown in sec. However, it may be 
used even beyond the adiabatic limit as the exact correlation function that extrapolates be- 
tween asymptotic regions, in which particles are well-defined. One interesting feature of the 
solitonic gauge fields is that the back-reaction from pairs through the induced electric field 
is smaller than the background gauge fields by order of the fine structure constant. There- 
fore, the back-reaction problem can be avoided for the solitonic gauge fields. It remains 
an open question to show the existence of solitonic gauge field including the back-reaction, 
for which the quantum master equation for the correlation functions is now coupled to the 
induced electric field, which in turn coupled to the correlation functions themselves. 
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